Generating and stabilizing the GHZ state in circuit QED: Joint measurement, Zeno effect and 

feedback 
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In solid-state circuit QED system, we extend the previous study of generating and stabilizing two-qubit Bell 
state [Phys. Rev. A 82, 032335 (2010)], to three-qubit GHZ state. In dispersive regime, we employ the 
homodyne joint readout for multiple qubits to infer the state for further processing, and in particular use it 
to stabilize the state directly by means of an alternate-flip-interrupted Zeno (AFIZ) scheme. Moreover, the 
state-of-the-art feedback action based on the filtered current enables not only a deterministic generation of the 
pre-GHZ state in the initial stage, but also a fast recovery from occasional error in the later stabilization process. 
We show that the proposed scheme can maintain the state with high fidelity if the efficient quantum measurement 
and rapid single-qubit rotations are available. 

PACS numbers: 03.67.Bg 32.80.Qk 42.50.Lc 42.50.Pq 



Quantum entanglement is novel and useful because it exhibits 
correlations that have no classical analog, and is one of the 
key ingredients for quantum technology applications such as 
quantum teleportation, quantum cryptography, quantum dense 
coding, and quantum computation |l|]. Conventional approach 
to generate quantum entanglement is via unitary two-qubit 
gates, using the necessary qubit-qubit interactions. Interest- 
ingly, instead of employing entangling gates, one can exploit 
quantum measurement as an alternative means to achieve sim- 
ilar goal. In both cavity and circuit QED systems, promising 
ideas along this route were proposed to probabilistically cre- 
ate entangled states by means of the homodyne measurement 
alone Hi]. 

The circuit QED system iQ-Jst], a solid-state analog of 
the conventional quantum optics cavity QED, is a promising 
solid-state quantum computing architecture. This architecture 
couples superconducting electronic circuit elements, which 
serve as the qubits, to harmonic oscillator modes of a mi- 
crowave resonator, which serve as a "quantum bus" that medi- 
ates inter-qubit coupling and facilitates quantum measurement 
for the qubit state. Moreover, quantum measurement in this 
system can be carried out by operating in the dispersive limit, 
i.e., with a detuning between the resonator and the qubit much 
larger than their coupling strength. In this limit, the qubit- 
resonator interaction induces a qubit-state-dependent shift on 
the resonator's frequency. Then, by measuring the resonator 
output voltage with a homodyne scheme, information about 
the qubit state is obtained. 

With these advantages, the schemes proposed in Refs. iSHt] 
for using measurement to create, respectively, two- and three- 
qubit entangled states in the circuit QED system, are attrac- 
tive. However, in addition to the drawback of being proba- 
bilistic, the measurement-only approach cannot stabilize the 
generated state. To resolve this problem, the technique of 
quantum feedback control may emerge as a possible route 
In a recent work IIQ] , a feedback scheme was analyzed for the 
creation and stabilization of the two-qubit Bell states. Owing 
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to using the dispersive joint-readout of multiple qubits and 
performing proper feedback, the scheme leads to results su- 
perior to some previous ones. For instance, it enhances the 
concurrence to values higher than 0.9, by noting the 0.31 ob- 
tained in Ref. ifTHl : and also it avoids the experimental diffi- 
culty in the jump-based feedback [l3], or the complexity for 
a state-estimation feedback 10, [H, 14ll . 

In this work, we extend the study of generating and stabi- 
lizing the two-qubit Bell state in Ref. iflOll to the three-qubit 
GHZ state. To our know ledg e, unlike the feedback control of 
two-qubit Bell state ifiol - fiafiel inil, schemes for stabilizing 
three-qubit GHZ state are not well investigated. In Ref. ^ 
elegant analysis is carried out for preparing the Dicke state 
of an ensemble of atoms (qubits) in cavity, by means of pro- 
jective measurement in dispersive limit, and as well by us- 
ing feedback to make the scheme deterministic. However, the 
ability of stabilizing the generated Dicke state against deco- 
herence is not demonstrated. In the present work, we plan to 
exploit the advantages of the homodyne joirit readout in dis- 
persive regime for multiple qubits, to infer the state for further 
processing, and in particular to stabilize the target state di- 
rectly by means of an alternate-flip-interrupted Zeno (AFIZ) 
scheme. Also, the state-of-the-art feedback action properly 
designed according to the measurement current enables both 
a deterministic generation of the pre-GHZ state in the initial 
stage, and a fast recovery from occasional error in the later 
stabilization process. 

Before proceeding to the details of the proposed scheme, 
we first briefly outline the control efficiency. For the deter- 
ministic generation of the pre-GHZ state, keeping track of the 
joint measurement information together with simple individ- 
ual qubit rotations in our scheme will either lead to a direct 
subsequent success of target state generation with probabil- 
ity 1/2, being higher than 1/4 given by the naive rerunning 
scheme which additionally needs the difficult "data clearing" 
procedures, or avoid to clear the wrong state before rerunning 
the generation scheme. More importantly, the AFIZ stabiliza- 
tion protocol, based only on an alternate but regular qubit flips 
and continuous measurement, can maintain very high fidelity 
(higher than ninety percent) for considerably long time (much 
longer than the single-qubit decoherence time). In principle. 
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FIG. 1: Schematic diagram of the circuit QED with three qubits, 
together with a microwave transmission measurement in dispersive 
limit. The Cooper-pair box qubits are fabricated inside a supercon- 
ducting transmission-line resonator and are capacitively coupled to 
the voltage standing wave. 



if the continuous measurement can approach the effect of fast 
repeated strong projective measurement, i.e., the ideal Zeno 
effect, the AFIZ scheme can stabilize the pre-GHZ state for 
arbitrarily long time. Meanwhile, another remarkable advan- 
tage of the AFIZ scheme is that it greatly simplifies the unitary 
manipulations on qubits, compared to either the state-based 
or current-based continuous feedback. Moreover, associated 
with the AFIZ scheme, an auxiliary alarm to the failure of 
AFIZ control, which may occur occasionally owing to the fi- 
nite strength of the measurement, can restart the fast determin- 
istic generation of the target state by simply using two-qubit 
flips, one qubit {n/2) rotation, and two projective measure- 
ments. 

Model and Formalism. — We consider a specific circuit 
QED system consisting of three superconducting qubits cou- 
pled to the fundamental mode of a microwave resonator cav- 
ity. The qubits, the resonator cavity and their mutual cou- 
pling are well described by the Jaynes-Cummings Hamilto- 
nian SQ: 



H 



a),a^a + S(a^ + a) + ^ 



(1) 



The operators crT(cr+) and ' ) are, respectively, the lowering 
(raising) operators ^or the jth qubit and the resonator cavity 
photon (hereafter we would name it as cavity photon for sim- 
plicity), (jji- is the frequency of the cavity photon, and and 
gj are the jth qubit transition energy and its coupling strength 
to the cavity photon. In this work we consider a three-qubit 
setup as shown schematically in Fig. 1, which can result in the 
required sign, sgn{gi) - sgnigo) - -sgnigj). The £-term in 
Eq. ([T]) stands for a microwave driving to the resonator cav- 
ity that is employed for the task of measurement. Explicitly, 
S = eg-™™' + c.c, where the driving frequency can differ from 
the cavity photon frequency, i.e., s - w„, ^ 0. More- 
over, we will focus on a dispersive limit measurement 
i.e., with the energy detuning A^ - cOr - Q.j much larger than 
gj. In this limit, the canonical transformation, H^ff ^ U^HU, 
where U = exp[2j /lj(flcrt - a'crp] with Aj = gj/^j, yields 
(in a joint rotating frame with the driving frequency oj,,, with 



respect to both the cavity photon and each individual qubit) 

Heff ^ Arfl'q + [ea + ea') + y/t<jj +Xi)^ + X^^'"'"^^- 



(2) 



where ojj - Qj-ajm, Sindxj - s]l^j- Here we have neglected 
the virtual cavity-photon-mediated effective coupling between 
|[ubits, as is appropriate for sufficient detuning between qubits 

In the circuit-QED system, the measurement is typically 
performed via a homodyne detection of the transmitted mi- 
crowave photons. The photon's leakage from the resonator 
cavity is described by a Lindblad term K'D[a]p in the master 
equation, where k is the leakage rate and the Lindblad super- 
operator acting on the reduced density matrix p is defined by 
2D[a]p = flpa ' - j{a^'a,p]. However, conditioned on the out- 
put homodyne current, i.e., /hom(0 - K{a + a')c{t) + y/K^it), 
there will be an additional unravelling term in the conditional 
master equation, 'H[fl]pc^(f). Here, <(■ ■ ■ ))c(f) s Tr[(- ■ ■ )pcit)] 
with pc{t) the conditional density matrix, and 'H[a]pc = apc + 
Pea' - Tr[(fl + a')pc\pc- And, the Gaussian white noise ^(f), 
which has the ensemble average properties £[^(0] - and 
E[^{t)^(t')] - 6{t - f'), stems from the quantum-jump related 
stochastic nature. 

In this work, we assume a strongly damped resonator cav- 
ity, which enables to adiabatically eliminate the cavity photon 
degree of freedom SlUt]. Qualitatively, observing the effec- 
tive coupling Yj^j=\ Xp'^'^^j in Eq. dU, we can understand that 
the fluctuation of the photon number will cause a pure dephas- 
ing backaction onto the qubits, with a joint dephasing operator 
Jz = Y.]=i where 5j = Xjlx and^ = Ij%iXj/N. Thus, 
we can expect the following results after adiabatic elimination 
of the photon's degree of freedom: (i) - the dephasing term 
~ D[JJpc', ( ii) - the unravelling term ~ 'J-{[J,]pc^{t); and ( Hi) - 
the homodyne current /hoin(0 ~ {Jz}c{t) + ^{t)- Indeed, follow- 
ing the standard procedures of adiabatic elimination iH [l9ll . 
an effective quantum trajectory equation (QTE) involving only 
the degrees of freedom of qubits can be obtained as 



Pc = -Cpc + yD[J,]p, + ^9{[J,]pcm, (3) 
in which the Liouvillian is defined as 



■CPc 



j j 

Here we have also assumed a resonant driving, i.e.. A,- = 0. 
In Eq. (O, jj and y^j are the relaxation and dephasing rates 
caused by the surrounding environments. Since the external 
dephasing can be strongly suppressed by proper design of the 
superconducting qubits, we will thus neglect it in our follow- 
ing simulations. The ypj terms, with -/pj - kA^., stem from 
the so-called Purcell effect, describing an indirect qubit deco- 
herence induced by the damping of the cavity photons. Since 
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we assumed mutually distinct frequencies of qubits, we can 
treat the Purcell effect under secular approximation, i.e., ne- 
glecting the cross terms in 2)['EjgjO']]Pc- (The cross terms 
characterize interference between the radiation from differ- 
ent qubits, and become important only if the qubit frequen- 
cies are sufficiently close.) With these considerations, the 
Purcell-effect-induced and other environment-caused deco- 
herences can be equally treated. Therefore, we combine them 
by y = jj + jpj. Moreover, the measurement-backaction in- 
duced dephasing rate - 2i\a'^x^ Ik, with a - —2ielK. And 
finally, the information-gain rate Y,„ in Eq. ^ is in general re- 
lated to the backaction dephasing rate in terms of the quantum 
efficiency, 77 - Y,„l{2Vd). 

Filtering the Output Current. — After adiabatic elimi- 
nation of the cavity degree of freedom, we obtain an effec- 
tive measurement operator 7, - Yj^j=\ ^i^j which has a num- 
ber of discrete eigenvalues. However, in the practical homo- 
dyne measurement, the output record is the homodyne current 
which we rewrite as dihomit) - y/T^,{J;)c(t)dt + dW{t), where 
dW{t) - ^{t)dt is the Wiener increment that has the statistical 
properties E[dW{t)] = and E[dW{t)dW(s)] = 6(t - s)dt. In 
this context, we remind that we are actually employing a se- 
ries of photons by using their transmission through the cavity 
to measure the qubit state. Accordingly, the stochastic Wiener 
increment in the homodyne current just characterizes the very 
weak (partial) random collapse caused by the individual mea- 
suring photons. It is the sum of the {Jz)c{t) term plus the 
Wiener increment that corresponds to the measurement record 
after the (infinitesimal) partial collapse, while {Jz)c{t), from 
its definition {J-)c{t) - Tr:[J.pc{ty\, exposes the information of 
state before the infinitesimal collapse. 

Of great interest and very usefully, instead of knowing 
{Jz)c(t) from the usual ensemble measurement by repeating a 
large number of realizations (noting that the quantum mechan- 
ical expectation indicates a statistical or ensemble-average in- 
terpretation), we can approximately obtain it in the context 
of continuous weak measurement by averaging the homo- 
dyne current over a properly chosen window of time. This 
has some similarity to the ergodic assumption in statistical 
physics, where the average along time is assumed to replace 
the ensemble average. In practice, we can low-pass filter 
the homodyne current over a small time window ifTsll . and 

get a smoothed signal as 7hom(0 = jj J,_j e"''"*'"'''^f/hom(T), 
where is a low-pass filtering parameter, and the factor TV 
normalizes the smoothed signal to a maximum magnitude of 
unity. Numerical test shows that /hom(0 indeed coincides with 
{Jz)c{t) satisfactorily. Thus, in the numerical simulations of 
this work, we simply use {Jz)c{t) as a state indicator to guide 
our feedback manipulations. In particular, after the qubits 
have experienced sufficient transmission measurement by a 
large number of photons and fully collapse onto one eigen- 
state of the measurement operator, we can in principle unam- 
biguously infer it from the filtered current {Jz)c{t). 

Finally, we remark that, in regard to feedback design, 
{Jz)c(t) in the homodyne current dihomit) is informative about 
the state and is thus useful, while the dW{t) term is yet harm- 
ful. In other words, the dW{t) term in the homodyne current 



has distinct roles when using the current to update the state 
versus to perform feedback. In doing the former, it is neces- 
sary; while in doing the latter, it is useless and should better 
be erased. This understanding was demonstrated in the recent 
work ifioll . 

Deterministic Generation of the Pre-GHZ State. — Since 
our control target, say, the GHZ state |000) + |111), can be 
easily obtained by a simple flip of single qubit (i.e. the third 
one) from |001) + |110) which may be accordingly named 
as pre-GHZ state, we can aim to control this pre-GHZ state 
instead. From the QTE formalism outlined above in disper- 
sive regime, we see that an effective measurement operator 
reads 7- - Y^j ^i^y Noting that the consequence of a quan- 
tum measurement is to collapse an arbitrary state onto one of 
the eigenstates of the measurement operator, we may design 
7; = cr^j + cr^ + 2cr^, which makes the pre-GHZ state be its 
eigenstate with eigenvalue /. - 0. In practice, this can be re- 
alized by setting the dispersive shifts : xi '■ X'i - ^ '■ 1:2. 

More specifically, let us start with an initially separable 
state; 

|*P,) = ^(|0) + |l))i®^(|0> + |l))2®^(|0) + |l»3 

V2 V2 V2 

= 4=[ 1000) + 1111) + (1010) + 1100)) 

+ (1011) + 1101)) + (1001) + 1110))]. (4) 

Performing the above designed homodyne measurement, in 
an individual single realization, would collapse I*!*,) stochasti- 
cally onto one of the eigenstates of J.- According to the prin- 
ciple of quantum projective measurement, I*!*, ) would collapse 
onto the pre-GHZ state |001) + |110) with probability 1/4, as 
a result of getting record 7. - 0. However, there are probabil- 
ities to get other results. That is, the state would collapse onto 
1011) + 1101) or 1100) + 1010) with probability 1/4, depending 
on the result to be 7, = 2 or -2. It may also collapse onto 
|1 1 1) or |000) with probabiUty I /8 if one gets 7, = 4 or -4. 

What we described above is in fact a measurement alone 
scheme to generate the pre-GHZ state stochastically. Be- 
low we show that proper current-based feedback manipula- 
tions can make the scheme deterministic. If, unfortunately, 
we did not get the result 7, = 0, we may adopt the fol- 
lowing distinct procedures based on the specific measure- 
ment results obtained, (i) - If the result is 7; - 2, which 
indicates the state |011) + |101) projected out, we perform 
a (Ti-flip on the first qubit and a njl - cr,, rotation on the 
third qubit. Noting that |011) + |101) can be rewritten as 
(|01) + |10))® |1), it is clear that the above rotations will trans- 
form it to 1000) + 1111) + 1001) + 1110), which then has a 
new probability of 1 /2 in the successive measurement to be 
collapsed onto the pre-GHZ state |001) + |110). (ii) - Simi- 
larly, for the result 7j — -2, a cTt-flip on the first qubit and a 
37r/2-cr,, rotation on the third one can be performed to achieve 
the same goal as described in (i). ( Hi) - If the measurement re- 
sult is J, = 4 or -4, which indicates the state |111) or |000) 
obtained, we then apply a. njl — cTy or a 37r/2 - cr,, rotation 
on each qubit, making the state return back to the initial one 
(|0) + |l))i ® (|0) + |1))2 ® (|0) + |1))3 that allows to rerun the 
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FIG. 2: (color online) Two representative quantum trajectories show- 
ing the deterministic generation of the pre-GHZ state. 



generating procedures. 

In Fig. 2 we show two representative quantum trajectories 
which are deterministically guided to the pre-GHZ state. The 
quantity {J;)c{t) plotted here is an appropriate indicator for 
the finally collapsed state. Here, we remark that our present 
scheme of entangled state generation is efficient. Generally 
speaking, probabilistic scheme of quantum information pro- 
cessing needs to recycle the process with the same initial state. 
This will require a procedure of clearing the unwanted data 
(qubit states), which is, unfortunately, not an easy job in quan- 
tum system compared to its classical counterpart. On the con- 
trary, in addition to that we do not require "data clearing", 
our scheme in case (ii) only needs two single-bit rotations and 
the subsequent success probability of projection (i.e., 1/2) is 
higher than the recycling scheme which has a success prob- 
ability 1/4. Even in the worst case (Hi), the manipulation of 
qubit rotations is identical to rerunning the generation proce- 
dure, but the process of "data-clearing" is avoided. 

Stabilization Using Flip-Interrupted Zeno Projection 
and Feedback. — In the above, we discuss an efficient deter- 
ministic scheme to generate the pre-GHZ state. However, un- 
der the unavoidable influence of the surrounding environment, 
this state will degrade if we do not provide proper active pro- 
tections. Below we first propose an alternate-flip- interrupted 
Zeno (AFIZ) stabilization scheme to enhance considerably the 
lifetime of the pre-GHZ state, then design additional manip- 
ulation to prevent the state from degradation. Actually, the 
pre-GHZ state is an eigenstate of the measurement operator 
J^. Then, if one performs a continuous observation (measure- 
ment) on it, the quantum Zeno effect would attempt to freeze 
the state. Using Zeno effect to stabilize quantum state is an in- 
teresting topic in quantum physics and particularly in quantum 
computing ll20[|2lll . However, as we will see shortly, general- 
ization of the Zeno protection from single to multiple qubits 
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FIG. 3: (color online) (a) State fidelity under the conventional quan- 
tum Zeno (not the AFIZ) stabilization for the pre-GHZ state, showing 
the result (solid line) much better than the uncontrolled one (dashed 
line), (b) Detailed inspection for the Zeno pulled-back state in (a), 
|>P(f)) - a(t)\00l) + yS(f)|110), showing a gradual deviation from the 
target state \^t) = (|001) + |110))/ VI. (c) Unconscious output cur- 
rent for the changing state IVit)). Single qubit decoherence rate: 
7 = O.OIFrf. 



will suffer more complexities. As a major contribution of the 
present work, we will see that the proposed AFIZ stabilization 
scheme can greatly improve the control quality. 

In Fig. 3(a) we plot the state fidelity under Zeno protection 
against the one in the absence of such protection, where the 
Zeno effect is automatically realized via the continuous J-- 
type measurement as discussed above. We see that, provided 
the measurement strength is much stronger than the decoher- 
ence rate, the effect of Zeno stabilization is obvious. Regard- 
ing the underlying mechanism of Zeno stabilization for multi- 
ple qubit state, it is analogous to that for single qubit. That is, 
while the environment is causing the state to leave away from 
the target state to other unwanted states, the relatively strong 
continuous measurement is at the same time pulling it back. 
Since our target state is a superposition of |001) and |110), 
it should, however, be more fragile than the one-component 
state in regard to the Zeno stabilization. 

We may understand this point better as follows. Owing to 
coupling with environment, the qubits would experience an 
entangling evolution with the environment, for an infinites- 
imal interval of time. Then, the measurement projects the 
qubit state back. Since each individual qubit couples to the 
environment independently, this pull-back action via measure- 
ment from entangling with the environment cannot guaran- 
tee the reorganized components |001) and |110) with exactly 
unchanged superposition weights. And, unfortunately, the 
changed superposition amplitudes cannot be distinguished by 
the measurement, since the arbitrary superposition of |001) 
and 1 110) is the eigenstate of the measurement operator and 
will result in the same output current. In Fig. 3(b) and (c) we 
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show, respectively, the gradual change of the superposition 
amplitudes away from the initial value 1 / V2, and the corre- 
sponding output current. 

The Zeno effect induced pull-back action, in terms of 
quantum measurement language, corresponds to a null-result 
record of spontaneous emission of the qubits. Conditioned on 
the null result of spontaneous emission, an effective Hamilto- 
nian governing the qubits state evolution reads H^^ - H^a - 
'5 Zij=i ^^j^j^ where H^a stands for the qubits Hamiltonian 
contained in Eq. (|2]i. With this effective Hamiltonian acting 
on the target pre-GHZ state \^t) = (|001> + |110))/ V2, or 
more generally on I*?) = fflOOl) + jS|l 10), leads to an effective 
evolution as |>P(f)> = (ae-y'l^\QQl) + ySe-'"|110))/|| ■ ||, where 
II ■ II denotes the normalization factor. We have checked that 
this effective evolution is in perfect agreement with the re- 
sults from numerical simulation, e.g., the one shown in Fig. 
3(b). Then, favorably, we get an insight that the change of the 
superposition amplitudes is owing to the unbalance of qubit 
state "0" and "1" in the components |001) and |110). Based 
on this observation, quite simply, if we flip simultaneously 
each of all the three qubits after the Zeno stabilization contin- 
uing for a time period r, the pre-GHZ state will be restored 
from the flipped state after an equal time interval r. In prac- 
tice, the time interval t can be chosen as a few F^', since in 
such a short timescale the variations of the state amplitudes 
a and p are negligibly small, provided y « F^. We can ex- 
pect and will demonstrate in the following that this alternate 
evolution, which is named in this work as an alternate -flip- 
interrupted Zeno (AFIZ) scheme, is capable of stabilizing the 
pre-GHZ state very efficiently. As an interesting remark, this 
AFIZ scheme has certain similarity to the refocus technique in 
the echo physics in quantum optics or nuclear magnetic res- 
onance, where similar flip is manipulated to cause an inverse 
evolution towards the initial state, i.e., to produce an echo. 

Unfortunately, to implement the above AFIZ protection in 
practice, there should exist very small but nonzero probabili- 
ties to collapse the qubits state to |000) and a mixture of 1 100) 
and |010), owing to finite strength of the measurement, (i) 
- For the result of |000), the output current would trigger a 
feedback action as described in the deterministic generation 
scheme, which will force the state rapidly back to the pre- 
GHZ state under the guided efficient projection of measure- 
ment, (ii) - For the mixture of |100) and |010), the output 
current will also trigger a feedback action as described in the 
deterministic generation scheme. As a result, besides being 
projected to |000) and |111) which will be further guided to 
the pre-GHZ state, a mixed state with |001) and |110) will 
be filtered out by the measurement. Very unfortunately, this 
mixed state is not the pre-GHZ state, but with the same zero 
output current. To eliminate this error, one can perform a flip 
action on the third qubit. Then, a mixed state with |000) and 
1 1 1 1) is formed and the rapid deterministic generation proce- 
dures will be triggered. 

In Fig. 4 we show the numerical result of stabilizing the pre- 
GHZ state, based on the measurement and feedback schemes 
described above. The stabilization dynamics is illustrated by 
both the state fidelity and the output current. We see that at 
most times the current is zero, only interrupted occasionally 
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no. 4: (color online) (a) Fidelity of the pre-GHZ state under the 
combined AFIZ-plus-feedback stabilization, (b) The corresponding 
output current. Single qubit decoherence rate: y = O.OlFj. 
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FIG. 5: (color online) (a) Average fidelity of the pre-GHZ state over 
1000 quantum trajectories, (b) Fidelity of an individual realization 
with 7 = O.OOlFrf, showing perfect control result under this even 
weaker decoherence when compared to y = O.OIF,/ in Fig. 4. (c) The 
full state density matrix at a specific time in (b). 



by jumps between 0, +2 and +4. The flat current indicates 
the stage of Zeno stabilization, during which the quality of 
the pre-GHZ state is maintained at desirable high level, with 
a state fidelity larger than 0.9. 

The state fidelity of single quantum trajectory has certain 
stochasticity, owing to the measurement induced quantum 
jumps. We may follow the conventional way to employ the 
fidelity of the ensemble average state as a reliable figure of 
merit to characterize the control quality. Figure 5(a) shows 
the results of ensemble average fidelity, for y = lO^^F^ and 
lO^^F,/, respectively. We notice that, even for y = 10"~F^, 
which in most cases such as the one- or two-qubit feedback 
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control is taken as a tolerable error rate, the average fidelity 
can be higher than 0.9, while for smaller error rate such as 
J - lO^^Frf the average fidelity can reach nearly unity and 
the individual quantum trajectory also shows perfect control 
result as illustrated in Fig. 5(b). In addition, in Fig. 5(c) we 
present the full density matrix for a representative state during 
the Zeno stabilization stage, which clearly reveals the quality 
of the protected pre-GHZ state. 

Concluding Remarks. — Finally, we make a number of 
remarks before summarizing the work, {i) - Some approxi- 
mations are involved, for instance, the rotating-wave approxi- 
mation contained in the Jaynes-Cummings model Eq. (1), the 
effective Hamiltonian Eq. (2) in a dispersive regime, and the 
adiabatic elimination of cav ity p hotons leading to Eq. (3). In 
the recent work by Liu et al lUOtl . all these approximations are 
properly justified. ( ii) - The main problem of doing feedback 
in circuit QED is the lack of efficient homodyne detection. 
Currently, the way to perform homodyne and heterodyne de- 
tection is to first amplify the signal before mixing it on a non- 
linear circuit element of some kind. As a consequence, the 
extra noise added by the amplifier will reduce the quantum ef- 
ficiency and prohibit quantum limited feedback. It seems that 
this situation is to be changed quickly, for instance, by de- 
veloping Josephson parametric amplifiers which can be real- 
ized in superconducting circuits ll22ll . - In our numerical 
simulation, we did not explicitly include the nonunit quan- 
tum efficiency in the homodyne detection of the field. Af- 
ter adiabatic elimination of the cavity photon degree of free- 
dom, the nonunit quantum efficiency of homodyne detection 



will reduce the effective information-gain rate F,,, in Eq. (3). 
This implies an emergence of an extra non-unravelling de- 
phasing term in the quantum trajectory equation. However, 
for the present particular study, this term only results in de- 
phasing among the pre-GHZ state and the others in the single 
quantum-trajectory realizations. In addition to simple intu- 
itive expectation, we have numerically examined that lower- 
ing the quantum efficiency by some acceptable amount does 
not obviously change the results, (iv) - Experimental verifica- 
tion of the GHZ state is of great interest and is analyzed the- 
oretically in the recent work by Bishop et al IH. In order to 
observe a violation of the Bell-Mermin inequality, relatively 
high signal-to-noise ratio in performing the measurement is 
required, which is unfortunately beyond the existing scope of 
experiment. But, optimistically, the situation is expected to 
change in the near future by fast experimental progresses. 

To summarize, we have presented a promising quantum 
control scheme for deterministic generation and stabilization 
of three-qubit GHZ state in the solid-state circuit QED sys- 
tem. The scheme largely depends on a joint-readout of mul- 
tiple qubits in dispersive regime, which enables not only to 
infer the state for further processing, but also to stabilize the 
target state directly by means of an alternate-flip-interrupted 
Zeno (AFIZ) projection. The proposed scheme was demon- 
strated by quantum trajectory simulations which show satis- 
factory control effects. 
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